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ABSTRACT 



^ I We present 3.5m Apache Point Observator}o second-epoch spectra of four low- 

^ I metaUicity emission-hne dwarf galaxies discovered serendipitously in the Data 

|.^ I Release 5 of the Sloan Digital Sky Survey (SDSS) to have extraordinary large 

O ' broad Ha luminosities, ranging from 3 x 10^^ to 2 x 10^^ erg s~^. The oxygen 

Q I abundance in these galaxies is very low, varying in the range 12 + log 0/H = 7.36 

- 7.99. Such extraordinarily high broad Ha luminosities cannot be accounted for 
by massive stars at different stages of their evolution. By comparing with the 
^ ' first-epoch SDSS spectra, we find that the broad Ha luminosities have remained 

constant over a period of 3 - 7 years, which probably excludes type Iln supernovae 
as a possible mechanism of broad emission. The emission most likely comes from 
accretion disks around intermediate-mass black holes with lower mass limits in 
the range ~ 5x10^ Mq - 3x10^ Mq. If so, these four objects form a new class 
of very low-met allicity AGN that have been elusive until now. The absence of 
the strong high-ionization lines [Ne v] A3426 and He ll A4686 can be understood 
if the nonthermal radiation contributes less than ~ 10% of the total ionizing 
radiation. 



^The Apache Point Observatory 3.5-meter telescope is owned and operated by the Astrophysical Research 
Consortium. 
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Subject headings: galaxies: abundances — galaxies: irregular — galaxies: active 
- galaxies: ISM — H ll regions — ISM: kinematics and dynamics 



1. INTRODUCTION 

Active galactic nuclei (AGN) are thought to be powered by massive black holes at 
the centers of galaxies, accreting gas from their surroundings. Observations of AGN show 
that they generally possess a high metallicity, varyirig fro m solar to supersolar metallic- 



ities (jStorchi-Bergmann et al. I Il998l : iHamann et al. II2002I ). While the derived metallici 



ties do depend on the detailed model assumptions, this appears to be a solid conclusion. 
Gas metallicity is know n to be strongly correlated with the stellar mass of the host galaxy 
( JTremonti et al. 112004 ). Since AGN are usually found in massive, bulge-dominated galaxies 
that have converted most of their gas into stars by the present epoch, their gas metallicities 
are generally high. A question then arises: do low-meta llicity AGN ex i st? If so, can we find 
them in low-mass galaxies? To address these questions, iGroves et al. I (120061 ) have searched 
the Sloan Digital Sky Survey (SDSS) Data Release 4 (DR4) spectroscopic galaxy sample of 
over 500,000 objects to select out ~170,000 emission-line galaxies with high S/N spectra. 
They then use diagnostic line ratios to select out 23,000 Seyfert 2s galaxies. Imposing an 
upper mass limit of 10^° Mq to restrict themselves to low-mass galaxies, they are left with a 
sample of only ~ 40 AGN, which they found appear to have metallicities around half that of 
typical AGN, i.e. having solar or slightly s ubsolar values. The sa me high metallicity range 
is found in the samp le of low-mass AGN of I Greene fc Ho I (1200 7bl ). Assessing their findings. 
Groves et al. I (120061 ) are led to another question: "Why are there no AGN with even lower 
metallicities?" In this paper, we suggest that these low-metallicity AGN do exist although 
they are extremely rare. 

In the course of a long-range prograrn to search for extremely metal- deficient emission- 



line dwarf galaxies, llzotov et al. I (l2007bl ) have used the SDSS DR5 database of 675,000 
spectra to assemble a large sample of emission-line galaxies. Two criteria were applied: 
1) the [O III] A4363 line must be detected to allow for a direct determination of element 
abundances; and 2) obvi ous high-met all i city AGN spectra are exclude d. Thus, contrary to 
Groves et al. I (120061 ) and iGreene fc Ho I (l2007bl ). llzotov et al. I (1200 7bl ) were not specifically 
looking for AGN. These criteria resulted in a sample of ~10,000 emission-line galaxies (ELG). 
Whil e studying tha t sampl e to look for ELGs with broad components in their strong emission 
lines, llzotov et al. I (l2007bl ) came across four galaxies with very unusual spectra. The general 
characteristics of the four galaxies are given in Table [H Their absolute magnitudes are typical 
of dwarf galaxies. Because of their relatively large distance ( z ~ 0.1-0.3) and relatively small 
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angular sizes (~ l"-2", only slightly larger than the seeing disk), their SDSS images (Fig. 
[1]) do not show much details. They possess a compact structure. Two galaxies, J1025+1402 
and J1047+0739, have an approximately round shape, while the other two more distant 
galaxies, J0045+1339 and J1222+3602, have a distorted shape suggestive of mergers. Their 
colors are not blue like the other ELGs, but vary from red to yellow to green. Their spectra 
shown in Fig. [21 resemble those of moderately to very low-metallicity high-excitation H ll 
regions: their oxygen abundances are in the range 12+log 0/H ~7. 4-7.9, i.e. their heavy 
el ement mass frac ti ons v ary from Z(^/ 1Q to Z(^/5 if the sol ar calibration 12+logO/H=8.65 
of lAsplund et al. I (120051 ) is adopted. Ilzotov et al. I (l2007bl ) found that there is however a 
striking difference: the strong permitted emission lines, mainly the Ha A6863 line, show very 
prominent broad components. These are characterized by somewhat unusual properties: 1) 
their Ha full widths at zero intensity FWZI vary from 102 to 158 A, corresponding to 
expansion velocities between 2200 and 3500 km s~^; 2) the broad Ha luminosities L^^ are 
extraordinarily large, varying from 3 x10^^ to 2x10 * ^^ erg s ~^. This is to be compared with 
the range lO^^-lO^'' erg s~^ found by Ilzotov et al. I (J2007bl ) for the other ELGs with broad- 
line emission. The ratio of Ha flux in the broad component to that in the narrow component 
varies from 0.4 to 3.4, as compared to 0.01-0.4 for the other galaxies; 3) the Balmer lines 
show a very steep decrement, suggesting collisional excitation and that the broad emission 
comes from very dense gas (A'^e^lO^ cm~^). Evidently, these galaxies are exceedingly rare, 
since they constitute only 4/675,000 or 0.0006% of our original sample. 



To account for the broad line emission in these four objects, Ilzotov et al. I (l2007bl ) have 
considered various physical mechanisms: a) Wolf-Rayet (WR) stars; b) stellar winds from Ofp 
or luminous blue variable stars; c) single or multiple Supernova (SN) remnants propagating in 
the interstellar medium; d) SN bubbles; e) shocks propagating in the circumstellar envelopes 
of typeHn SNe; and f) AGN. While mechanisms a-d can account for L^r ~ 10^^ to lO^'' erg 
s~^, they cannot provide for luminosities that are 30 to 200 tim es greater. Th e se very large 
luminosities are more likely associated with SN shocks or AGN. Ilzotov et al. I (2007bt) have 
considered typelln SNe because their Ha luminosities are larger (~10^^-10^^ erg s~^) than 
those of the other SN types. Up and III, and they decrease less rapidly. 

To decide whether type Hn SNe or AGN are responsible for the broad emission in these 
galaxies, monitoring of their spectral features on the relatively long time scale of several 
years is necessary. If broad features are produced by Hn type SNe, then we would expect a 
decrease in the broad line luminosities. No significant temporal evolution would be expected 
in the case of an AGN. Additionally, higher signal-to-noise ratio spectra are necessary to 
put better constraints on the presence of the high-ionization [Ne v] A3426 and He ll A4686 
emission lines, good indicators of a source of hard non-thermal radiation. In order to check for 
temporal evolution, we have obtained second-epoch spectra of the above four galaxies with 
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broad emission, using the 3.5ni Apache Point Observatory (APO) telescopcl- We describe 
the observations in §2. We discuss in §3 the main properties of the broad emission and show 
that they can be accounted for by low-metaUicity intermediate-mass AGN. Our conclusions 
are summarized in S4. 



2. OBSERVATIONS 

New high signal-to-noise ratio optical spectra were obtained for the four galaxies listed 
in Table [H, using the 3.5 m APO telescope on the nights of 2007 November 15 and 2008 
February 6. The observations were made with the Dual Imaging Spectrograph (DIS) in the 
both the blue and red wavelength ranges. A l'.'5x360" slit was used. In the blue range, 
we use the B400 grating with a linear dispersion of 1.83 A/pix and a central wavelength of 
4400A, while in the red range we use the R300 grating with a linear dispersion 2.31 A/pix 
and a central wavelength of 7500A. The above instrumental set-up gave a spatial scale along 
the slit of 0'.'8 pixel"^ on the night of 2007 November 15 and of 0'.'4 pixel"^ on the night of 
2008 February 6 (the latter scale is smaller by a factor of 2 than the previous one because we 
did not perform pixel binning for the latter spectra), a spectral range ~3600 - 9600A and 
a spectral resolution of 7A (FWHM). The slit was oriented along the parallactic angle and 
the total exposure time was 45 minutes for each galaxy. The observations were broken up 
into 3 subexposures to allow for removal of cosmic rays. The Kitt Peak IRS spectroscopic 
standard stars Feige 110 (2007 November 15), Feige 34 and G191B2B (2008 February 6) 
were observed for flux calibration. Spectra of He-Ne-Ar comparison arcs were obtained at 
the beginning of each night for wavelength calibration. 



The data reduction procedures are the same as described in lThuan fc Izotovi (120051 ). The 
two-dimensional spectra were bias subtracted and fiat-field corrected using IRAF. We then 
use the IRAF software routines IDENTIFY, REIDENTIFY, FITCOORD, TRANSFORM 
to perform wavelength calibration and correct for distortion and tilt for each frame. Night 



^In fact. Ilzotov et al. I ( 2007a ) have identified 5 objects as AGN candidates. We have not included here 



the fifth candidate, J2230-0006=PHL 293B, because its broad Ka luminosity is only S.GxfO"^^ erg s ^, some 
10'^ - 10* times lower than those of the o ther four AGN candidates. The broad Ha luminosity of J2230-0006 



has been erroneously given in Table 8 of Ilzotov et al. I (J2007bl ) as having 10 times its true value. Given this 



relatively low luminosity and the fact that a new 3.5m APO spectrum of J2230-0006 shows that its broad 
hydrogen lines have a P Cygni profile with a blue-shifted absorption, the broad emission probably originates 
from a stellar wind rather than from an accretion disk around a AGN. Most likely, the broad emission in 
J2230-0006 is caused by a strong outburst in a bright luminous b l ue va riable (LBV) star. Similar broad 



hydrogen emission has been detected recently by iPustilnik et al. I (|2008| ) in the extremely metal-deficient 
dwarf galaxy DDO 68. 
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sky subtraction was performed using the routine BACKGROUND. The level of night sky 
emission was determined from the closest regions to the galaxy that are free of galaxian 
stellar and nebular line emission, as well as of emission from foreground and background 
sources. A one-dimensional spectrum was then extracted from the two-dimensional frame 
using the APALL routine. An extraction aperture of I'.'S x 4" was adopted. Its area is 
similar to that of the 3" round aperture used in the SDSS spectra and allows to gather > 
90% of the light of o ur compact galaxies. Comparison of the line intensities of the SDSS 
(jlzotov et al. I l2007bl ) and 3.5m spectra (Table [2]) shows that they agree well, validating 
our choice of the extraction aperture. Before extraction, the three distinct two-dimensional 
spectra of each object were carefully aligned using the spatial locations of the brightest part 
in each spectrum, so that spectra were extracted at the same positions in all subexposures. 
We then summed the individual spectra from each subexposure after manual removal of the 
cosmic rays hits. 

The resulting APO spectra of the four galaxies are shown in Fig. [2l A strong broad Ha 
emission line is present in all spectra, very si milar to the one see n in the SDSS spectra of 
the same galaxies obtained 3-7 years earlier (jlzotov et al. Il2007bl ). The broad components 
in the H/3 emission line are considerably weaker, suggesting a steep Balmer decrement and 
hence that the broad emission originates in a very dense gas. In the cas e of J1047+0739, 
broad He I is also present, as in the SDSS spectrum (jlzotov et al. Il2007bl ). 



3. RESULTS AND DISCUSSION 



3.1. Element abundances 



We have derived element abundances from the narrow emission line fluxes. These fluxes 
have been measured using Gaussian fitting with the IRAF SPLO T routine. They have 
been corrected for both extinction, using the reddening curve of (jWhitford I Il958l ). and 
underlying h ydrogen stellar abso rption, derived simultaneously by an iterative procedure as 
described by llzotov et al. I (jl994j ) and using the observed decrements of the narrow hydrogen 
Balmer B.S A4101, H7 A4340, H/3 A4861 and Ha A6563 lines. It is assumed in this procedure 
that hydrogen line emission is produced only by spontaneous transitions in recombination 
cascades, i.e. we neglect possible collisional excitation. Such a situation usuall y holds in low 



density H ll regions ionized by stellar radiation such as those considered here. Ilzotov et al. 



(J2007al ) have shown that the use of different reddening curves has little influence on the 
extinction-corrected fluxes (relative to the H/5 flux) in the optical range and modify only 
slightly the extinction coefficient. The reason is that the spectra are corrected in such a 
way so that the relative intensities of the extinction-corrected hydrogen lines correspond 
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to their theoretical values. The extinction-corrected fluxes 100 x /(A) //(H/?) of the narrow 
lines for each galaxy, together with the extinction coefficient C(H/5), the equivalent width 
of the H/3 emission line EW(H/?), the H/? observed flux F(H/5) and the equivalent widths 
of the underlying hydrogen absorption lines EW(abs) are given in Table [21 The physical 
conditions and element abundances of the H ll regions in the four AGN candidates are 
derived from the narrow line fluxes fo llowing llzotov et al. I ( 20061 ). The element abundances 
derived from the 2.5m SDSS spectra (llzotov et al. Il2007bl ) and from the 3.5 m spectra (this 
paper) are shown in Table [31 We find the abundances and the abundance ratios obtained 
for the four galaxies from the two sets of observations agree well. The derive d abundances 



are in the range characteristic of low-metallicity dwarf emission-line galaxies (llzotov et al 



20061 ). This implies that the narrow emission lines arise primarily from regions ionized by 
stellar ionizing radiation, and that any contribution of a possible AGN component to the 
narrow line emission is small. 



3.2. Broad emission and diagnostic diagrams 



Using Gaussian fitting, we have measured the fiuxes of the broad component of the 
Ha line after subtraction of the narrow component. These fiuxes are shown in Table HI 
Comparison of SDSS and APO broad Ha fiuxes shows that they have remained nearly 
constant (with variations < 20%) over a period of ~ 3 - 7 years. This likely rules out 
the hypothesis that the broad line fiuxes are due to type Iln SN because their Ha fiuxes 
should have decreased significantly over this time interval. T here is one known exception, 
SN 1996cr, where the Ha fiux has remained constant for 7 yrs (JBauer et al. II2008I ). although 
here the luminosity of Ha is only ~ 10^^ erg s~^. Maintaining an Ha luminosity of 10^^ - 
Xg42 gj,g g-i ^Qj, such durations in our objects would require nearly the entire energy budget 
of a SN (> 10^° - 10^^ ergs). We thus rule out type Iln SNe. 

There remains the AGN scenario. Can accretion disks around black holes in these low- 
metallicity dwarf galaxies account for their spectral properties? The spectra of the four 
objects do not show clear evidence for the presence of an intense source of hard nonthermal 
radiation: the [Ne v] A3426, [O ll] A3727, He ll A4686, [O l] A6300, [N ll] A6583, and [S 
ll] AA6717, 6731 emission lines, which are usually found in the spectra of AGN, are weak 
or not detected. Aside from He ll A4686, the apparent weakness of such emission lines, 
however, may be accounted for by the low metallicities of our galaxies. None of our objects 
were detected in the NVSS or FIRST radio catalogs, demonstrating that they are faint radio 
sources. Another way to check for the presence o f an AGN in a galaxy is to check for its 



location in the emission-line diagnostic diagram of [Baldwin et al. I (Il98l[ ) (BPT). It can be 
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seen that all four objects lie in the region corresponding to star- forming galaxies (SFG), to 
the left of the region occupied by AG N with low-mass black holes and with metallicities 
ranging from 2 to 1/4 that of the Sun (JGreene &: Ho Il2007bl ). However, their locations in 
the SFG region do not necessarily disqualify them as AGN candidates. Photoionization 
models of AGN show that lowering their meta llicity moves them to the left of the BPT 
diagram, so that they end up in the SFG region (JGroves et al. II2006I : IStasihska et al. II2006I ). 
Thus the BPT diagram is unable to distinguish between SFGs and low-met allicity AGN. 
Admitting that there is an AGN in our dwarf galaxies, can we account for the weakness 
of the high-ionization lines? Photoionization models with only AGN nonthermal ionizing 
radiation do predict detectable He ll A4686 and [Ne v] A3426 emission lines. To make the 
observed spectra agree with the models, one solution is to dilute the nonthermal ionizing 
radiation from the AGN by thermal radiation from s urrounding hot massive stars. In F ig. 
E^, we show the results of our CLOUDY calculations (lFerlandlll996l : iFerland et al.lll998l ) of 
H II regions ionized by a composite radiation consisting of different proportions of stellar 
and nonthermal radiation. Two curves, characterized by different metallicities, are shown by 
solid lines: the lower one is for 12+logO/H=7.3 and the upper one is for 12+logO/H=7.8, 
typical of the metallicities of our objects. Each model point is labeled by the ratio R of 
nonthermal-to-thermal ionizing radiation. A slope a = -1 has been adopted for the non- 
thermal power-law spectrum over the whole wavelength range under consideration {fi, oc z/"). 
The calculations have been done with a number of ionizing photons (5t/i=10^^ s~^ for stellar 
radiation, Qnonth = RQth for nonthermal radiation and N(.=1Q^ cm~^. Higher densities would 
move the curves to the righ t. For the ionizing stellar radiation, we adopt Costar models by 
Schaerer fc de Koter I (119971 ) with a heavy element mass fraction Z = 0.004 and an effective 



temperature of 53,000K corresponding to a starburst age of < 3 Myr. The dotted hues in Fig. 
|3^ show the corresponding models with a = -2. They are very similar to the models with a 
= -1 when R < 1, but fall below for R> 1. It is seen that models with 12+logO/H=7.8 and 
in which the nonthermal ionizing radiation contributes ^10% of the total ionizing radiation 
can account well for the location of all four galaxies in the BPT diagram, independently of 
the slope of the power-law spectrum. 

How about the high-ionization lines? Can their absence be due to high dust extinction 
in the central part of the galaxy? We consider this possibility unlikely because in this case, 
broad Ha emission from the accretion disk surrounding the black hole would not be seen 
also. We turn next to the properties of the ionizing spectrum for an explanation. In Fig. [3]d, 
we show the diagnostic diagram for [Ne v] A3426/H/3 vs. [N ll] A6583/Hq; (thick lines) and 
He II A4686/H/3 vs. [N ll] A6583/Ha (thin lines). As in Fig. Oi, CLOUDY models with a = 
-1 and -2 are shown by solid and dotted lines. They are also characterized by (5t/i=10 



53 



for stellar radiation, Q 



nonth 



RQth for nonthermal radiation and A^e=10^ cm ^. Higher 



densities would shift curves to the right. The vertical dashed line separates models with 
12+logO/H = 7.3 (Fig. [3b, left) from those with 12+logO/H = 7.8 (Fig. [3)d, right). The 
shaded rectangle shows the region of the upper limits of ~ 1% - 2% of the B.f3 flux, set for 
[Ne v] A3426/H/3 and He ll A4686/H/3 in our objects, in the observed range of their [N ll] 
A6583/Ha ratio. There are several points to be made concerning Fig. [3]d. First, in contrast 
to the low- ionization line fluxes (Fig. [3^), the predicted fluxes of the high-ionization lines 
are more strongly dependent on the slope of the power-law spectrum of the nonthermal 
radiation. Second, while there is no significant change of the predicted fluxes of the He ll 
A4686 emission line with metallicity, the predicted flux of the [Ne v] A3426 linearly scales 
with 12-|-logO/H. Third, at a given metallicity, the detectability of the high-ionization lines 
depends on two parameters: one is the ratio R of nonthermal-to-thermal ionizing radiation 
and the other is the slope of the power-law ionizing spectrum. If we adopt 12+logO/H = 7.8 
as typical for our galaxies, then Fig. [3)d shows that models that satisfy the non-detectability 
limit of the high-ionization lines (i.e. that fall within the shaded box) are characterized 
either by a steep slope and a not excessively small R {a = -2 and i? ~ 0.1) or by a shallower 
slope and a very low R {a = -1 and R ~ 0.03). We expect intermediate-mass black holes 
to have accretion disks that are hotter than those around supermassive black holes, and 
hence their ionizing spectrum to have a shallower slope, nearer -1 rather than -2. If that is 
the case, then the fraction of ionizing nonthermal to thermal radiation is very small in our 
dwarf galaxies, ~ 3%. It is also possible that the absence of strong high-ionization lines is 
caused by a high covering factor of the accretion disk. In this case the hard radiation would 
be absorbed inside the dense accretion disk and no high-ionization forbidden lines would be 
formed. 



3.3. Black hole virial masses 

We now estimate the masses of the central black holes. It has been shown (see e.g. 



Kaspi et al. 1 12000| ) that continuum and broad line luminosities in AGN can be used to 



determine the size and geometry of the broad emission-line re gion and the mass o f the 



central black hole. Examining a large sample of broad-line AGN, I Greene fc Ho I (120051 ) have 
found that the Ha luminosity scales almost linearly with the optical continuum luminosity 
and that a strong correlation exists between the Ha and H/3 line widths. On the basis of 
these two empirical correlations, those authors have derived the following relations for the 
central black hole mass: 
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M, 



BH 



4.4 X 10 



L 



5100 



0.64 



/FWHMh/3 



M 



0) 



(2) 



where Lhq and L^ioq are respectively the broad Ha and continuum (at A = 5100A) lumi- 
nosities, and FWHMh/3 and FWHMHa are respectively the full widths at half maximum of 
the H/? and Ha emission lines. 

In Table E] we list the extinction-corrected broad Ha luminosities L(Ha) and continuum 
luminosities ALa(5100) for the four galaxies, as derived from the SDSS spectra. The extinc- 
tion coefficient was set equal to the one derived for the narrow Balmer hydrogen lines. Since 
the reddening due to dust extinction in dense regions may be larger than that derived from 
the narrow hydrogen emission lines, the derived L(Ha) should be considered as lower limits. 
The L(Ha) and ALa(5100) of ou r galaxies follow clos ely the correlation between Ha and 
continuum luminosities found by I Greene &: Ho I (120051 ). This implies th a t our galaxies are 
very likely the same type of objects as those considered by I Greene fc Ho I (120051 ). Therefore, 
we can use Eqs. [1] - [2] for the determination of the central black hole masses. For the flux 
and luminosity determinations, we first fit the line profiles by a single Gaussian. However, 
we find that single Gaussian fits do not reproduce well the broad low-intensity wings of the 
Ha line. This suggests that that the line broadening is caused not only by gas motions but 
also by light scattering. This hypothesis is supported by a steep Balmer decrement (the 
Ha-to-H/9 fiux ratio is greater than ~ 7 in our four galaxies), implying high gas densities 
and possibly high optical depths in the Ha line. Therefore, in addition to a single Gaussian 
profile fit, we have also fitted the Ha line by a Voigt profile which is a superposition of both a 
Gaussian profile and a Lorentzian profile. In Table [5], we show the FWHMs of the Gaussian 
components in both cases. The FWHMs are smaller by a factor of < 2 in the case of the 
Voigt profile fitting as compared to the single Gaussian fitting. For the central black hole 
mass determination, we have adopted the Gaussian FWHMs derived from the more accu- 
rate Voigt profile fit ting. We have also set FWHMh/j to be equal to FWHMho, following 



Greene fc Ho I (120051 ). The masses MBH(Ha) and Mbh(5100) derived from the broad Ha 
and continuum luminosities are shown in Table O The masses derived from the two meth- 
ods are in good agreement, as expected from the relation between L(Ha) and ALa(5100) 
for our object s and the correlation between the Ha and continuum luminosities found by 



Greene fc Ho I (120051 ). The derived masses of the central black holes in our galaxies are in 
the range ~ 5x1 0^ M(^ - 3x10^ Mp,, lo wer or similar to the mean black hole mass of 1.3 x 
10^ Mq found by iGreene fc Ho I (1200 7bl ) for their sample of low- mass black holes. Since the 



luminosities used to derive the masses of the central black holes are lower limits, the derived 
masses should also be considered as lower limits. 
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4. CONCLUSIONS 

We study here the broad hne emission in four low-metalhcity star-forming dwarf galaxies 
with 12+logO/H ~ 7.4 - 8.0, i.e. with metalhcities between 1/19 and 1/5 that of the Sun. 
We have arrived at the following conclusions: 

1. The steep Balmer decrements of the broad hydrogen lines and the very high lumi- 
nosities of the broad Ha line in all four galaxies (3x10^^ to 2x10^^ erg s~^) suggest that the 
broad emission arises from very dense and high luminosity regions such as those associated 
with supernovae (SNe) of type Iln or with accretion disks around black holes. However, the 
relative constancy of the broad Ha luminosities over a period of 3-7 years likely rules out the 
SN mechanism. Thus, the emission of broad hydrogen lines is most likely associated with ac- 
cretion disks around black holes. If so, these four objects would harbor a new class of AGN 
that are extremely rare (in 0.0006% of all galaxies). These AGN would be intermediate- 
mass black holes residing in low-metallicity dwarf galaxies, with an oxygen abundance that 
is considerably lower than the solar or super-solar metallicity of a typical AGN. 

2. There is no obvious spectroscopic evidence for the presence of a source of a non- 
thermal hard ionizing radiation in all four galaxies: high-ionization emission lines such as 
He II A4686 and [Ne v] A3426 emission lines were not detected at the level < 1 - 2 percent 
of the H/3 flux. We have calculated a series of CLOUDY models with ionizing spectra which 
include both thermal stellar and nonthermal power-law ionizing radiation in order to account 
for the absence of the high-ionization lines. We find that the predicted fluxes of the high- 
ionization lines are below the detectability level if the spectral energy distribution /j, oc z^" 
of the ionizing nonthermal radiation has a ~ -1 and the nonthermal ionizing radiation is 
significantly diluted by the thermal stellar ionizing radiation contributing < 3 percent of the 
total ionizing radiation, or the ionizing spectrum is steeper {a ~ -2), and the nonthermal 
ionizing radiation contributes < 10 percent of the total ionizing radiation. 

3. The lower limits of the masses of the central black holes Mbh of ~ 5x10^ Mq - 
3x10^ Mq in our galaxies are among the lowest found thus far for AGN. 
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J0045+1339 for us. We thank Franz Bauer and John Hawley for useful discussions. Luis Ho 
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Table 1. General characteristics of low-metallicity AGN candidates 



Object R.A. (J2000.0) DEC. (J2000.0) Redshift g Mg 



SDSSJ0045+1339 


00 45 29.2 


+13 39 09 


0.29522 


21.80 


-18.56 


SDSSJ1025+1402 


10 25 30.3 


+14 02 07 


0.10067 


20.36 


-17.66 


SDSSJ1047+0739 


10 47 55.9 


+07 39 51 


0.16828 


19.91 


-19.23 


SDSSJ1222+3602 


12 22 45.7 


+36 02 18 


0.30112 


21.30 


-19.10 
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Table 2. Intensities of narrow emission lines from 3.5m spectra 



Ion 


100x/(A)/J(H/3) 


100x/(A)/J(H/3) 


100x/(A)/J(H/3) 


100x/(A)//(H/3) 






Galaxy 






J0045+1339 


J1025+1402 


J1047+0739 


J1222±3602 


3727 [0 II 


65.45 ± 2.47 


70.06 ± 6.30 


94.18 ± 2.20 


41.36 ± 2.27 


3868 [Ne in 


50.42 ± 2.05 


56.01 ± 5.29 


52.64 ± 1.35 


64.17 ± 2.68 


3889 He I + H8 




20.40 ± 4.21 


18.17 ± 1.02 


12.74 ± 3.01 


3968 [Ne ill] + H7 


29.75 ± 2.51 


40.16 ± 4.88 


33.31 ± 1.18 


47.17 ± 3.50 


4026 He I 


1.18 ±0.65 








4101 H(5 


25.04 ± 2.84 


24.47 ± 4.42 


26.78 ± 1.10 


26.95 ± 2.80 


4340 H7 


43.49 ± 2.49 


50.40 ± 5.80 


45.90 ± 1.30 


44.38 ± 4.01 


4363 [0 III 


14.78 ± 1.59 


24.44 ± 2.98 


10.68 ± 0.50 


20.77 ± 1.76 


4471 He I 






5.03 ± 0.33 




4861 H/5 


100.00 ± 3.12 


100.00 ± 13.09 


100.00 ± 2.24 


100.00 ± 6.33 


4959 [0 III 


244.61 ± 5.94 


200.48 ± 14.85 


183.36 ± 3.48 


283.29 ± 7.73 


5007 [0 III 


747.92 ± 17.18 


575.78 ± 34.53 


555.65 ± 10.10 


853.85 ± 22.21 


5876 He I 






20.38 ± 0.62 




6563 Ha 


279.78 ± 13.73 


271.10 ± 82.71 


281.42 ± 9.24 


277.54 ± 14.82 


6583 [N II 


6.54 ± 1.03 


4.06 ± 1.18 


4.50 ± 0.09 




6678 He I 






6.15 ±0.27 




6717 [S II 






6.53 ± 0.26 




6731 [S II 






4.23 ± 0.24 




7065 He I 






28.89 ± 0.74 




7136 [Ar III 






5.00 ± 0.26 




C{W) 


0.135 


0.110 


0.400 


0.100 


EW(H/3) A 


394 


45 


116 


196 


F{B.13Y 


0.12 


0.16 


0.70 


0.14 


EW(abs) A 


5.55 


0.10 


5.35 


5.50 



'^in units 10 ^^ erg s ^ cm ^. 
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Table 3. Element abundances 



Object 2.5m SDSS 3.5m APO 



12+log 0/H log N/0 log Ne/0 12+log 0/H log N/0 log Ne/0 

J0045+1339 7.94 ± 0.08 -1.23 ± 0.17 -0.77 ± 0.15 7.86 ± 0.05 -1.14 ± 0.08 -0.83 ± 0.07 

J1025+1402 7.36 ± 0.08 -1.16 ± 0.30 -0.86 ± 0.14 7.48 ± 0.08 -1.30 ± 0.14 -0.74 ± 0.12 

J1047+0739 7.99 ± 0.05 -1.21 ± 0.09 -0.73 ± 0.09 7.85 ± 0.02 -1.46 ± 0.03 -0.68 ± 0.04 

J1222+3602 7.94 ± 0.09 -1.02 ± 0.24 -0.84 ± 0.14 7.88 ± 0.05 • • • -0.70 ± 0.07 
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Table 4. Fluxes of the broad Ha emission line at different epochs 







2.5m SDSS 


3.5m APO 


Object 


T a 
-'br 


hr/Inar^ Date of Obs. 


T b 
-'br 


Date of Obs. 


J0045+1339 


16.4±1.7 


41.65 2000 Jan 12 


18.0±0.7 


2007 Nov 15 


J1025+1402 


165.0±5.7 


337.49 2004 Mar 11 


192.5±0.8 


2008 Feb 06 


J1047+0739 


289.2±9.1 


122.34 2003 Jan 31 


224.2±2.1 


2008 Feb 06 


J1222+3602 


16.1±1.8 


71.69 2005 Mar 13 


22.4±1.3 


2008 Feb 06 



^In units 10 ^^ erg s ^ cm ^. Flux errors are derived taking into account photon 
statistics in non-flux calibrated spectra. 



b 



In percent. 
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Table 5. Ha and continuum luminosities and masses of the black lioles^ 



Object Lbr(H.a)^ FWHMfer(Ha) ALa(5100)^ MBH(Ha)^ MbhISIOO)^ 



J0045+1339 2.74x10^^ 

J1025+1402 3.21x10^^ 

J1047+0739 1.57x10^2 

J1222+3602 2.80x10^^ 



Gaussian^' 


Voigt^'*^ 








1700 


1540 


7.59x10^2 


2.43x10^ 


2.00x10*^ 


1580 


680 


3.62x10*^2 


5.07x10^ 


2.50x10^ 


1920 


1050 


2.32x10^^ 


3.05x10^ 


1.91x10*' 


1600 


790 


7.17x10^2 


6.34x10^ 


5.10x10^ 



^Parameters are derived from the 2.5 m SDSS spectra (jlzotov et al. Il2007bl ). 

'^In units erg s~^. 

'^In units km s~^. 

^Derived by fitting a Gaussian profile. 

''Derived by fitting a Voigt profile, a combination of Gaussian and Lorentzian profiles. 
The FWHM is that of the Gaussian profile. 

^In solar masses. 
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J0045+1339 



»i^» km. 



m 



J102 5+1402 




J 1047+0739 




J 1222+3602 



Fig. 1. — 13" X 13" SDSS images of low-metallicity AGN. The angular diameters of the objects 
vary between 1" and 2", barely larger than the seeing disk. 
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Fig. 2. — Redshift-corrected 3.5 m Apache Observatory second-epoch spectra of four low- 
metalhcy emission-hne dwarf galaxies thought to contain AGN. The locations of the non- 
detected [Ne v] A3426 and He ll A4686 high-ionization emission lines are shown in all panels. 
Other emission lines are labeled in panel c). 
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Fig. 3. — a) The BPT diagram flBaldwin et al. 
ted are the ~1 0,000 ELGs fro m 



Izotov et aL 



19811) for low-ionization emission hnes. Plot- 



(l2007bl ) (cloud of points), the low-mass black 



hole sample of I Greene fc Ho I (l2007bl ) (asterisks) and the four low-metallicity AGN in Table 
[1] (stars') ■ The dashed line separ ates star-forming galaxies (SFG) from active galactic nuclei 
(AGN) (JKauffmann et al. II2003I ). The solid and dotted lines connect CLOUDY photoiniza- 
tion models computed for H ll regions ionized by a composite radiation consisting of different 
proportions of stellar and nonthermal radiation. The two upper curves are characterized by 
12+log 0/H = 7.8 and the two lower ones by 12+log 0/H = 7.3. Each model point is labeled 
by the ratio of nonthermal-to-thermal ionizing radiation. All curves have been calculated 
adopting a number of ionizing photons Q = 10^^ s~^ for stellar radiation, different slopes of 
the nonthermal spectral energy distributions fi, oc z^" (solid lines are for a = -1 and dotted 
lines are for a = -2). A density N^. = 10^ cm~^ is adopted. Higher densities would shift 
the curves to the right, b) The diagnostic diagram for high- ionization emission lines: [Ne v] 
A3426/H/3 vs. [N ll] A6583/Ha (thick lines) and He ll A4686/H/? vs. [N ll] A6583/Ha (thin 
lines). The same CLOUDY models as in a) are shown. The shaded region shows the upper 
intensity limits of the high-ionization lines [Ne v] A3426 and He ll A4686 relative to H/3 in the 
four galaxies considered here. The dashed vertical line separates models with 12+logO/H = 
7.30 (left) from those with 12+logO/H = 7.80 (right). 



